Anomalous Paramagnetic Effects in the Mixed State of LuNi 2 B 2 C 
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Anomalous paramagnetic effects in dc magnetization were observed in the mixed state of 
LuNi2B2C, unlike any reported previously. It appears as a kink-like feature for H > 30 kOe 
and becomes more prominent with increasing field. A specific heat anomaly at the corresponding 
temperature suggests that the magnetization anomaly is due to a true bulk transition. A magnetic 
flux transition from a square to an hexagonal lattice is consistent with the anomaly. 



In cuprate (high-T c ) superconductors, high-transition 
temperatures (T c ) and short coherence lengths (£) lead 
to large thermal fluctuation effects, opening a possibil- 
ity for melting of the flux line lattice (FLL) at temper- 
atures well below the superconducting transition tem- 
perature. A discontinuous step in dc magnetization and 
a sudden, kink-like drop in resistivity signified the first 
order nature of the melting transition from the vortex 
lattice into a liquidi^ In conventional type II supercon- 
ductors, with modest transition temperatures and large 
coherence lengths, vortex melting is also expected to oc- 
cur in a very limited part of the phase diagram^ but it 
has yet to be observed experimentally. In the rare-earth 
nickel borocarbides RNi 2 B 2 C (R = Y, Dy, Ho, Er, Tm, 
Lu), the coherence lengths (£ = 10 2 A) and superconduct- 
ing transition temperatures (16.1 K for R = Lu) lie be- 
tween these extremes, suggesting that the vortex melting 
will be observable and may provide further information 
on vortex dynamics. Indeed, Mun et al£ reported the 
observation of vortex melting in YNI2B2C, based on a 
sharp, kink-like drop in electrical resistivity. 

Recently, a magnetic field-driven FLL transition has 
been observed in the tetragonal borocarbides &1&*L The 
transition from square to hexagonal vortex lattice occurs 
due to the competition between sources of anisotropy and 
vortex- vortex interactions. The repulsive nature of the 
vortex interaction favors the hexagonal Abrikosov lat- 
tice, whose vortex spacing is larger than that of a square 
lattice. The competing anisotropy, which favors a square 
lattice, can be due to lattice effects (fourfold Fermi sur- 
face anisotropy) ^fi unconventional superconducting or- 
der parameter or an interplay of the twoi^iA In com- 
bination with non-negligible fluctuation effects, the com- 
petition leads to unique vortex dynamics right below the 
H c 2 line in the borocarbides, namely a reentrant vor- 



tex lattice transition^ Fluctuation effects near the upper 
critical field line wash out the anisotropy effect, stabiliz- 
ing the Abrikosov hexagonal lattice li^i^ Here, we report 
the first observation of paramagnetic effects in the dc 
magnetization M of the mixed state of LUN12B2C. The 
kink-like feature in M and the corresponding specific heat 
feature for H > 30 kOe signify the reentrant FLL tran- 
sition, which is consistent with the low-field FLL tran- 
sition line inferred from small angle neutron scattering 
(SANS) .2 

Single crystals of LUM2B2C were grown in a M2B flux 
as described elsewhere^ and were post-growth annealed 
at T = 1000 C° for 100 hours under high vacuum, typi- 
cally low 10~ 6 Torriii Samples subjected to a preparation 
process such as grinding, were annealed again at the same 
condition as the post-growth annealing. A Quantum De- 
sign magnetic property measurement system (MPMS) 
was used to measure ac and dc magnetization while the 
heat capacity option of a Quantum Design physical prop- 
erty measurement system (PPMS) was used for specific 
heat measurements. Electrical resistivity was measured 
by using a Linear Research ac resistance bridge (LR-700) 
in combination with a PPMS. 

The detailed dc magnetization of LUNI2B2C reveals an 
anomalous paramagnetic effect for H >30 kOe, where 
the magnetic response deviates from a monotonic de- 
crease and starts to rise, showing decreased diamag- 
netic response. The in-phase and out-of-phase compo- 
nents of the ac susceptibility Xac show a dip and the 
specific heat data show a feature at the corresponding 
temperature, reminiscent of vortex melting in high-T c 
cuprates^ Electrical transport measurements, however, 
do not exhibit any feature corresponding to the param- 
agnetic effect; e.g., a sharp drop in the electrical resis- 
tivity. The zero-resistance transition, rather, occurs at a 
much higher temperature, suggesting that the paramag- 
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FIG. 1: (Color online). Top panel: dc magnetization M(T) 
as a function of temperature at 20, 25, 30, 35, 40, and 
45 kOe. Dashed lines are fits from the standard local London 
model— and solid lines from the non-local London model— 
(see text). Bottom panel: temperature derivation of M(T) at 
corresponding magnetic fields. Arrows indicate the points 
where kink-like features start to appear 



netic effect is not related to vortex melting. It is instead 
consistent with a topological FLL change between square 
and hexagonal structures. 

The top panel of Fig. 1 shows dc magnetization M 
as a function of temperature at several fields. For H > 
30 kOe, kink-like features appear, which are marked by 
arrows. The anomalous increase can be easily seen as a 
sharp drop in dM/dT (arrows in the bottom panel). The 
magnetization reported here is independent of time and 
has no hysteresis between zero-field cool (ZFC) and field 
cool (FC) data within experimental accuracy, indicating 
that the measured value is an equilibrium magnetization. 
In the top panel of Fig. 1, we compare the data and some 
model calculations. Dashed lines are predictions from the 
standard local London model: 19 



4ttM = M ln(r)H c2 /H) 



(1) 



where Mo = (/)q/8ttX 2 , 77 is a constant of order unity, A 
the penetration depth, and (f>o the flux quantum. In the 
fit, H c2 was determined from our resistivity data (see 
Fig. 2) and M from H c2 with n = A/£ = 15. In or- 



der to get the best result, the fitting parameter -q was 
varied between 0.95 and 0.97 and the absolute ampli- 
tude of Mo was changed as a function of magnetic field. 
The local London model explains the monotonic decrease 
with decreasing temperature, but the fit becomes worse 
at higher field. In a clean system like LuNi2B2C where 
the electronic mean free path is long compared to the 
coherence length £o> the current at a point depends on 
magnetic fields within a characteristic length p, or nonlo- 
cal radius. Taking into account the nonlocal current-field 
relation in superconductors, a non-local London model 
was suggested— i 

- 4ttM = M [ln(l + H Q /H) - H Q /(H Q + H) + A], (2) 

where H a = 4> / (Air 2 p 2 ) and A = rn - ln(l + H a /r) 2 H c2 ) 
with 771 and 772 being order of unity. It is worth not- 
ing that the scaling parameter H c2 in the local theory is 
replaced by Hq in the non-local model. The nonlocal ra- 
dius p slowly decreases with increasing temperature and 
is suppressed strongly by scattering. The solid lines are 
best results from the model calculation where we used the 
temperature dependence of Hq and A from the literature 
for YNi2B2Ci2i Both the local and the non-local models 
explain the temperature dependence of M at low fields, 
while only the non-local model can describe the data at 
and above 35 kOe. The good fit from the nonlocal model 
is consistent with the equilibrium magnetization analy- 
sis of YM2B2C, 21 suggesting the importance of nonlocal 
effects in the magnetization. The many fitting param- 
eters in both fits, however, prevent us from making a 
definite conclusion as to which model better describes 
M(T). Nevertheless, we can extract the important con- 
clusion that the kink-like feature in the mixed state is a 
new phenomenon that needs further explanation. 

In the early stage of high-T c cuprate research, anoma- 
lous paramagnetic effects in M(T) were reported in the 
irreversible region and this effect was later attributed 
to the field inhomogeneity of the measured scan length 
in a SQUID magnetometer^ We tested various scan 
lengths from 1.8 cm to 6 cm for which the field inho- 
mogeneity varies from 0.005 % to 1.4 % along the scan 
length and found negligible dependence on the measur- 
ing length, which suggests that field inhomogeneity is 
not the source of the anomaly. A more definitive test 
used a conventional type II superconductor NbSe2 in a 
similar configuration. There was no such anomalies in 
NbSe2 as in the borocarbide. Taken together, we con- 
clude that the reversible paramagnetic effects are intrin- 
sic to LUN12B2C. We also emphasize that the the phe- 
nomena is different from the paramagnetic Meissner ef- 
fect (PME) or Wohlleben effect where the FC x be- 
comes positive whereas the ZFC x remains negative. The 
PME is an irreversible effect and occurs in the Meissner 
state, while the subject of this study is a reversible effect 
and takes place in the mixed state. 

In the top panel of Fig. 2 the reversible magnetiza- 
tion M (left axis) and the out-of-phase component of ac 
susceptibility x'ac (right axis) are shown as a function of 
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FIG. 2: (Color online). Top panel: Magnetization (left axis) 
and the imaginary part of the ac susceptibility (right axis) at 
40 kOe. Bottom panel: Resistive superconducting transition 
at 40 kOe (circles) and 50 kOe (squares). The lines are guide 
to eyes. 



temperature at 40 kOe. A dip appears both in x'ac ( n °t 
shown) and in x'ac a ^ the same temperature where M 
shows the paramagnetic anomaly. Since a dip in \ac is 
often related to vortex melting, it is natural to consider 
the vortex phase change from liquid to lattice or glass 
as a possible explanation. The resistive superconducting 
transition at 40 kOe (circles) and 50 kOe (squares) are 
shown in the bottom panel of Fig. 2. A resistive slope 
change in the transition region, that can be considered as 
a signature of the vortex melting^ was observed at 8.1 K 
and 6.6 K for 40 kOe and 50 kOe, respectively. The R = 
transition temperature, however, is much higher than the 
temperature where the dip occurs in Xac which argues 
against the vortex melting scenario as the physical origin 
of the anomalous paramagnetic effects. The increase in 
M at the transition temperature is also opposite from 
the decrease in the vortex melting interpretation. 

Recently, a structural phase transition in the FLL was 
suggested to explain another peak effect observed be- 
low the vortex melting line in YBCOi22i2£ The vanish- 
ing of a squash elastic mode gives rise to a topolog- 
ical FLL transition and leads to the new peak effect, 
while the softening of the shear modes cqq is relevant 
to the conventional peak effect in high-T c cuprates. 26 
The observation of the dip effect well below the melting 
line in LuNi2B>2C indicates that the anomalous param- 
agnetic effects are related to a change in the FLL and 
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FIG. 3: (Color online). Temperature dependence of A(4nM) 
at 30, 35, 40, and 45 kOe, where AM = M - M (nonlocal) . 
Inset: The peak values of A(4nM) are compared with those 
estimated from Eq. (3). 



the increase in M is also consistent with the FLL change 
where the Abrikosov geometrical factor (3 changes ,A2i2£ 
Fig. 3 shows the temperature dependence of the para- 
magnetic anomaly, A(4irM), at several magnetic fields, 
where A(4-7rM) is the magnetization after subtracting 
the monotonic, diamagnetic background obtained from 
Eq. (2). With increasing field, the peak becomes en- 
hanced and an additional peak is observed at 45 kOe. In 
extreme type II material (k >> 1/V2), the magnetiza- 
tion change due to a FLL transition is written as 



A(4ttM) = 



H C 2 — H 
2k 2 - 1 



1 

Pa 



1 

Jo 



(3) 



where (3a ~ 1.16 for a hexagonal FLL and /3<> w 1.18 for 
a square FLLpi^ In the inset of Fig. 3, we compared the 
peak intensity of A(47rM) (left axis) and the estimation 
from Eq. 3 (right axis) with k = 15. It is encourag- 
ing to see that the simple model qualitatively reproduces 
the field dependence of the paramagnetic contribution. 
However, the quantitative difference in absolute values 
suggests that a more elaborate model is required. 

The H — T phase diagram is shown in Fig. 4. The 
upper critical field line H C 2 was determined from the 
R = superconducting transition and is consistent with 
the temperature where Xac starts to have a non-zero 
value. The Hi line in the mixed state is the point where 
the reversible magnetization shows the paramagnetic ef- 
fects and H2 is the 2nd anomaly that appears above 
45 kOe (see Fig. 1). According to the Gurevich-Kogan 
non-local London model^ the anisotropic nonlocal po- 
tential, which is responsible for the low-field FLL transi- 
tion observed in SANSjS is averaged out by thermal vor- 
tex fluctuations near H C 2 . Since the interaction becomes 
isotropic, the hexagonal Abrikosov lattice is preferable, 
leading to the second FLL transition from square back 
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FIG. 4: (Color online). H — T phase diagram. The Hi line 
indicates the temperature where M(T) deviates from a mono- 
tonic decrease and the H2 line is where the second anomaly 
occurs above 45 kOe. The dashed line depicts qualitatively 
what fluctuation models based on the nonlocal London— or 
extended GL— ■ predict when the nonlocal radius is approxi- 
mately 2 — 3 £o- The dotted line is from theories without fluc- 
tuation effects?— *— The diamonds are SANS data, where the 
numbers next to the data indicate the degrees of the azimuthal 
splitting with which the transition line is determined. 9 The 
square and rhombic shapes are forms of vortex lattices. 



to rhombic (triangular) lattice as the field gets closer to 
the H C 2 line. The reentrant transition is predicted to 
occur well below the vortex melting line because the am- 
plitude of vortex fluctuations required to wash out the 
nonlocal effects is much smaller than that for the vortex 
melting. This prediction is consistent with our observa- 
tion that the Hi line is much below the R = tran- 
sition line. The dashed line depicts qualitatively what 
the fluctuations models based on the nonlocal London 14 
or extended GL— i predict, which nicely explains the Hi 
line. The dotted line is the FLL transition line that meets 
the H C 2 line both in the non-local London models and 
in the extended Ginzburg-Landau (GL) theory^ with- 
out fluctuation effects. We note that a direct compar- 
ison between the SANS and our data is difficult even 
though they are qualitatively similar. Since Hi heavily 
depends on the sample purity,™ a factor of 2 or more dif- 
ference in Hi has been easily observed even among pure 
compoundsi^iiS Further, the Hi line from the SANS also 
depends on the criteria used for the FLL transition (see 
Fig. 4). 

For vortex melting, where the lattice changes to a liq- 
uid, the transition involves latent heat and the specific 
heat shows a sharp peak at the transition temperature— ^ 
For a structural change in the vortex lattice, the transi- 
tion is probably of 2nd order because an infinitesimally 
small change of the angle (3 between adjacent vortex lines 
changes the symmetry. Based on the paramagnetic jump 



FIG. 5: (Color online). Specific heat difference C(H)/T - 
C{9T)/T vs T at 45 kOe for H || [001]. Inset: blow-up of the 
main panel at around 6.1 K. Different symbols correspond to 
different sets of measurements and attest to the reproducibil- 
ity of these results. 



(Fig. 1) and dHi/dT < (Fig. 3), Ehrenfest's relation 
at constant field predicts a suppression of C/T as the 
FLL changes from a rhombic to a square lattice. Fig. 5 
shows the specific heat data of LuNi2B2C at 45 kOe as 
a function of temperature. In addition to the supercon- 
ducting transition between 6.22 and 6.89 K, an anomaly 
is, indeed, observed at 6.12 K which corresponds to the 
anomalous paramagnetic effects. Depending on the back 
ground we choose, however, the anomaly can be consid- 
ered as either a suppression or a jumpjS Similar features 
at 40 and 50 kOe were also observed at the temperatures 
corresponding to the paramagnetic effects in H\. More 
sensitive measurements such as ac calorimctry will help 
in resolving the issue. Finally, we note that we are not 
able to discern any corresponding feature to the H2 line 
in C p or in Xac- More work is in progress to understand 
the second paramagnetic jump in M which appears for 
H > 45 kOe. 

In summary, we report the first observation of an 
anomalous paramagnetic jump in the magnetization of 
the mixed state of LUN12B2C. A dip appears in Xac at the 
same temperature as the paramagnetic effects, suggesting 
the relevance of the flux line lattice. The H — T phase di- 
agram is consistent with a FLL structural transition from 
square to hexagonal lattice just below the upper critical 
field line. The observation of an additional feature in 
the specific heat data at the corresponding temperature 
underscores the interpretation of paramagnetic effects as 
due to a reentrant FLL transition in LuNi2B2C. 
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